Electric Dipole Induced Spin Resonance in Quantum Dots 
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An alternating electric field, applied to a quantum dot, couples to the electron spin via the spin- 
orbit interaction. We analyze different types of spin-orbit coupling known in the literature and find 
two efficient mechanisms of spin control in quantum dots. The linear in momentum Dresselhaus 
and Rashba spin-orbit couplings give rise to a fully transverse effective magnetic field in the pres- 
ence of a Zeeman splitting at lowest order in the spin-orbit interaction. The cubic in momentum 
Dresselhaus terms are efficient in a quantum dot with non-harmonic confining potential and give 
rise to a spin-electric coupling proportional to the orbital magnetic field. We derive an effective spin 
Hamiltonian, which can be used to implement spin manipulation on a timescale of 10 ns with the 
current experimental setups. 

PACS numbers: 73.23.Hk, 73.63.Kv 



CZ2 



-a 
c 
o 

(N 
> 

o 

\o 
o 

-I— > 
c3 



i 

C 

O 

o 



X 



I. INTRODUCTION 

Coherent manipulation of electron spin is at the heart 
of spintronics 1 ^ and quantum computing with spins£ In 
the proposal of Ref. y, the spin of an electron confined 
to a quantum dot is used as qubit to store and process 
quantum information. A quantum register consisting of 
an array of such spin- 1/2 quantum dots is operated by a 
set of quantum gates that act on single spins and pairs of 
neighboring spins^ Among the simplest quantum gates 
are the spin rotations on the Bloch sphere. With the 
help of only a static magnetic field and an electron-spin- 
resonance (ESR) pulse, one can change the state of the 
spin qubit at will. It is important, however, that the ESR 
pulse can be applied locally to each of the quantum dots, 
ensuring that the spins are accessed independently from 
one another. For an ESR 4-5 to occur, usually, the elec- 
tron is exposed to an alternating magnetic field of a fre- 
quency uj ac that matches the electron Zeeman splitting. 
However, because strong local electric fields are easier to 
obtain than strong local magnetic fields, interest arises 
in spin resonance induced by electric fields. 

Recently, Kato et al&> have demonstrated three- 
dimensional control of spins in a GaAs/ALjGai_ x As het- 
erostructure with the use of an alternating electric field. 
The mechanism of spin coupling to the electric field re- 
lies on a specially engineered Lande g tensor in the het- 
erostructure, achieved by modulating the Al content dur- 
ing the MBE growthii The resulting g tensor is both 
anisotropic and space-dependent, and allows control over 
the direction and magnitude of the spin precession fre- 
quencyiS^ A g-factor modulation resonance (g-TMR) oc- 
curs similarly to an ESR, when the frequency of the 
electric field matches the Zeeman splitting. 6 Rashba and 
Efros 9 have further proposed to use the standard (Dres- 
selhaus 10 and Rashbaiii) spin-orbit couplings to achieve 
an electric dipole induced spin resonance (EDSR) in 
quantum wells. Rashba and Efros 9 have shown that 
the EDSR is highly efficient in quantum wells, promis- 
ing electron spin control on a timescale w^ 1 ~ 100 ps, 
where lor is the Rabi frequency. 4 These results have im- 



portant practical implications in spintronics, where spins 
of extended electrons are used as a resource to accom- 
plish information processing. In the context of quantum 
computing, however, interest arises in the spin of a lo- 
calized electron. A natural question is, therefore, "What 
is the microscopic mechanism of EDSR in quantum dots 
and how strong is the EDSR effect?" 

EDSR has nearly half-a-century long history. It has 
been first observed for extended electrons in bulk semi- 
conductors^ 4 ^ and studied more recently for donor- 
bound electrons in Cdi-^Mn^Se 1 ^ and extended elec- 
trons in two-dimensional electron gases 1 ^ 1 ^ and epilay- 
ersiii The "forbidden" electric-dipole transition between 
the electron spin-up and spin-down states becomes possi- 
ble in the presence of spin-orbit interaction. Absorption 
spectra of EDSR provide information about the value of 
the electron g factor and the strength of the spin-orbit 
coupling. In two-dimensional electron systems, one ex- 
pects the Dresselhaus spin-orbit interactioni - to be en- 
hanced compared to bulk semiconductors, because of the 
confinement of electron motion in one direction. Further- 
more, the Rashba spin-orbit interactional arises in het- 
erostructures lacking inversion symmetry, such as, e.g., 
heterojunctions. In some systems, the Rashba coupling 
constant can be efficiently tuned by electric fields^ 

In quantum dots^ the spin-orbit interaction is gener- 
ally suppressed due to complete localization of electron 
motion* 20 Typically, the quantum dot lateral size \<i is 
smaller than the spin-orbit length Xso, and any effect 
of the spin-orbit interaction is suppressed as a power of 
Ad/Aso and therefore is expected to be weak. This ex- 
pectation contrasts with the case of electrons in quantum 
wells, where the EDSR meets most favorable conditions^ 
The Zeeman interaction in quantum dots plays an impor- 
tant role for observing spin- orbit effectsiSii 2 ^ Without the 
Zeeman interaction, the Rashba and linear in momentum 
Dresselhaus spin-orbit terms do not contribute to spin- 
related phenomena at the first order of spin-orbit interac- 
tion. This "absence of spin-orbit" at the leading order in 
quantum dots has been discussed extensively in the liter- 
ature. 20 ! 21 : 22 Below, we show that a similar result arises 
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FIG. 1: Schematic of a setup for electric field control of spin 
via the spin-orbit interaction. The quantum dot (QD) con- 
tains a single electron with spin S = (h/2)a, deep in the 
Coulomb blockade valley. The gates 1 and 2 are used to gen- 
erate an alternating electric field E(t), which acts via the 
spin-orbit interaction on the electron spin. As a result, an 
electric dipole spin resonance (EDSR) occurs if the frequency 
of E(t) is tuned to match the Larmor frequency u>z = Ez/h. 



also for the cubic in momentum Dresselhaus terms in the 
case when the dot confining potential is quadratic and 
the perturbation is linear in the electron coordinates. 

In this paper, we consider the use of EDSR for control 
of individual electron spins in quantum dots. We derive 
an effective spin Hamiltonian for a quantum-dot electron, 
subject to ac electric fields. We show that there are two 
major mechanisms of EDSR in quantum dots. One arises 
from the linear in momentum Dresselhaus and Rashba 
spin-orbit couplings in combination with the Zeeman in- 
teraction. The other arises from the cubic Dresselhaus 
terms in combination with the cyclotron frequency. We 
estimate the strengths of both EDSR effects and compare 
them to the ordinary ESR. We find that despite a strong 
suppression, compared to quantum wells, the EDSR in 
quantum dots is still an efficient mechanism of spin ma- 
nipulation and can be used alone or together with ESR 
to achieve control of spin on a timescale iv^ 1 ~ 10ns. 



II. EDSR SETUP 

We consider a quantum dot containing a single elec- 
tron with charge — e and spin S = (K/2)(T, where cr = 
cr y , <t z ) are the Pauli matrices. The quantum dot is in 
the Coulomb blockade regime with extraction (E/_) and 
addition (J7+) energies that are large compared to the 
temperature, so that the dot occupation remains con- 
stant. An external electric field E(r,t) is applied to the 
quantum dot. In practice, E(r,t) can be generated by a 
pair of gates, as sketched in Fig. ^ to which an ac signal 
of frequency uj ac is supplied from an external circuit (not 
shown). The Hamiltonian describing the quantum dot 
electron in the external alternating field reads 



H = H + V(r,t), 



(1) 



where V(r, t) — e f r dr' ■ E(r' , t) is the potential energy 
of the electron in the external electric field and Hq is the 



"unperturbed" Hamiltonian (see further). In particular, 
for an electric field constant in space E(r, t) = E(t), the 
potential energy reads V(r 7 t) — eE(t) ■ r. 

For practical applications, it is a good idea to use two 
gates, as shown in Fig. ^ because this allows larger am- 
plitudes of E(t) to be applied to the quantum dot, while 
still maintaining the dot within the same Coulomb block- 
ade valley. Ideally, the ac voltage drop is distributed 
between the two gates symmetrically and the dot po- 
tential is kept constant by counteracting potential shifts 
quadratic in the electric field. For a harmonic quantum 
dot, the desired ac potential reads 



V(r,t) = eE(t) ■ r 



(eE(t)) 2 
2m e u) 2 



(2) 



where m e is the electron effective mass and u>o is the 
oscillator frequency. Then the only effect of the ac signal 
on the dot confinement is shifting the dot center as a 
function of time by the amount 



ro(t) 



eE{t) 



(3) 



The amplitude of »"o(i) is going to be a relevant parame- 
ter in our following analysis. Therefore, setups in which 
the dot can be easily moved on the substrate by gates 
are particularly interesting in the context of this paper. 
We discuss the case of ro ~ Xso in Sec. IVII whereas for 
the bulk of the paper we restrict ourself to r <C Xso ■ 
The Hamiltonian Hq consists of several terms, 



Hq — Hd + Hz + H, 



so, 

where Hd is the Hamiltonian of a confined electron, 

,,2 



H d 



P 
2m e 



U(r), 



(4) 



(5) 



with p = —ihd/dr + (e/c)A(r) being the electron mo- 
mentum, c the speed of light in vacuum, and U(r) the 
quantum dot confining potential. We restrict our consid- 
eration to quantum dots with strong confinement along 
one axis, such as, e.g., quantum dots defined in a two- 
dimensional electron gas (2DEG). For GaAs, the 2DEG 
lies, typically, in the crystallographic plane (001) and has 
a width d ~ 5 nm, which ensures a strong size quantiza- 
tion along z || [001]. The in-plane motion of the electron 
is described by the Hamiltonian (JSJ, where r = (x,y) 
is the electron in-plane coordinate; whereas the trans- 
verse motion (along z) has already been integrated out in 
Eqs. CJ-®- In the absence of external time-dependent 
fields, A{r) accounts for the orbital effect of a static mag- 
netic field B. Assuming that B is constant in space, we 
have A(r) — B z (—y/2,x/2,0) in the symmetric gauge. 
Note that the in-plane components B x and B y are not 
present in A(r), because the motion along z is strongly 
quantized {d <§; y j / hc/eB x ^) . 

The magnetic field B also induces a Zeeman splitting 
Ez = 9^-bB and a spin quantization axis n = B/B via 
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the Zeeman interaction, 



1 1 
Hz = 2 S ^ bB ' a = o zU ' <T ' 



(6) 



where g is the electron g-factor and /is is the Bohr mag- 
neton. In GaAs, the magnitude of the g- factor is anoma- 
lously small (g « —0.44) compared to other AmBy 
semiconductors. The Zeeman energy is, therefore, much 
smaller than the cyclotron energy Hui c — heB z /m e c by a 
factor gm e /m <C 1 (with m being the electron mass in 
vacuum), for magnetic fields applied transversely to the 
2DEG. In Sec.[V] we derive an efficient spin-electric cou- 
pling that is proportional to fkv c , but present only in non- 
harmonic quantum dots. We remark that the magnetic 
field is an important ingredient in our EDSR scheme, 
since at B = no spin-electric coupling can be obtained 
at the first order of the spin-orbit interaction (see fur- 
ther) . 

In Eq. (J3J, Hso stands for the spin-orbit Hamiltonian. 
We start with considering the so-called "linear in p" spin- 
orbit interaction, 



relativistic nature of B(t). Furthermore, the proxim- 
ity of the top gates to the 2DEG decrease the displace- 
ment current enclosed by the magnetic field lines pen- 
etrating the quantum dot. Using Eq. we estimate 
B = /jnc~ 1 L z dE/dt w 1CP 6 G, for a quantum dot which 
is L z = 100 nm below the gates plane and an electric field 
E(t) — Eosm(ui ac t), with amplitude (ikEq = 10 2 V/cm 
and frequency oj ac /2ir — 10 9 Hz. As we show below, a 
much stronger effective magnetic field arises from the 
EDSR effect in the present setup, and therefore, the dis- 
placement current can be safely ignored. 

Recently, a sizable ESR effect has been obtained with 
the help of a wire placed on top of a GaAs double dot. 27 
In this case, J 7^ and the magnitude of B(t) is estimated 
from Eq. © to be B — iKiic~ x Ij{L y + L z ), where I is 
current in the wire and 2L y is the lateral size of the wire. 
For L y = L z = 100 nm and 1=1 mA, the magnetic field 
obtained in this setup is on the order of B ~ 10 G. 



III. SPIN-ELECTRIC COUPLING 



Hso = a{p x a y - PyO x ) + l3(-p x cr x + p v cr v ), (7) 

which is the sum of the Rashba (a)!^ and 2D Dresselhaus 
(o '^iLiii spin-orbit interactions. This type of spin-orbit 
interaction gives rise to a sizable phonon-induced spin 
relaxation rate l/Tn^is— of the same order of magnitude 
as experimentally measured ^SS^i Moreover, in the 2D 
limit the linear in p spin-orbit interaction is dominant, 
because cx 1/d 2 . 

In the standard EDR scheme, an alternating magnetic 
field is generated by a current in a nearby conductor. In 
our setup (see Fig. , no charge flow is ideally present 
between the gates. However, the alternating electric field 
E(r,z,t) gives rise to a displacement current, with the 
current density 



J D {r,z,t) 



k dE(r, z, t) 



4tt 



dt 



(8) 



where k is the electric permittivity. The external mag- 
netic field B acquires, thus, an ac-component, B — > 
B + B(t), where B{t) = V x A(t). The vector potential 
A.(r, z, t) is obtained as usual from Ampere's law2£ 



V 2 A- 



(J + Jl 



(9) 



where fi is the magnetic permittivity and J is the charge 
flow density (in our case J = 0). In Eq. ©, we adopted 
the Coulomb gauge V • A. = and used the notation 
V = (d/dr,d/dz). 

The magnetic field B(t) couples to the electron spin via 
the Zeeman interaction in Eq. © [with B — * B + B(t)], 
giving rise to an ESR source, which can be used, in 
principle, for spin manipulation in quantum dots. How- 
ever, the amplitude of Bit) is, in practice, extremely 
small; it is proportional to 1/c, as expected from the 



Now, we focus on the electric-field component of 
the ac-signal and show that, together with the spin- 
orbit interaction Hso and Zeeman splitting Hz, it suf- 
fices to generate a sizable EDSR field in the quan- 
tum dot. For simplicity, we set Bit) — > from now 
on and choose A(t) = 0. As a result, we retain 
only a constant in space and time magnetic field B = 
B(cos ip sin sin ip sin -d, cos ?9) and an alternating electric 
field E(r,t) = (1/e) VU(r, t). Thus, we consider further 
the Hamiltonian in Eq. Q), assuming that Hq is time- 
independent and describes the dot in the absence of ac 
fields. 

We aim at diagonalizing Hq using a Schrieffer- Wolff 
transformation, similar to Ref. 0. We thus look for 
a transformation matrix S such that the transformed 
Hamiltonian Hq = cxp(S)Ho exp(— S) is fully diagonal, 
see Appendix lAl At B = 0, the ground state of Hq (and 
also of £?o) is a Kramers doublet, because the spin-orbit 
interaction is always time-reversal symmetric (at B = 0). 
We therefore choose to encode the qubit into the ground 
state Kramers doublet of the quantum dot. Owing to the 
mixed spin and orbital nature of the states an alternat- 
ing potential V(r,t), such as in Eq. (Q, couples to the 
qubit. We proceed to derive this coupling, by calculat- 
ing the transformation matrix S at the leading order of 
spin-orbit interaction, 



1 



L d + L 2 



-H 



so 



0{H, 



(10) 



where Ld and Lz are Liouville superoperators, i.e. 
L d A = [Ha, A] and L Z A = [H Z ,A], VA. The projec- 
tor V projects onto the diagonal (or degenerate) part of 
the Hilbert space of H& + Hz, which ensures applicabil- 
ity of "non-degenerate" perturbation theory. The cou- 
pling of spin to electric fields is then found by applying 
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the same Schrieffer- Wolff transformation to the potential 
V(r,t). We obtain the following effective Hamiltonian 
for our qubit in the presence of an alternating potential 
V(r, £), to leading order in the spin-orbit interaction, 



H, 



H z + (il>o\[S,V(r,t)]\rl>o), 



(11) 



where S is the transformation matrix in Eq. (|10() and 
\ipo) is the quantum dot ground state. For a quantum dot 
with a harmonic confining potential, U(r) = m P uj 'kr 2 / 2, 
the transformation matrix S was calculated in Ref. [23 to 
all orders of the Zeeman interaction and the first order of 
the "linear inp" spin-orbit interaction 10. For simplicity, 
we consider here only the linear in B terms, 



S = i£ 



m e u> 



n x 



C] 



(12) 



£ = (AlV.A^V.O), 

C = (\z 1 a/d 1 /,x+ 1 a/dx',o), 

where X± = h/m e ({3± a) are the spin-orbit lengths, and 
the vectors £ and £ are given in the coordinate frame x' = 
(x + y)/ \/2, y' = —(x — y)/ V2, and z' = z (see Fig. . 
The first term in Eq. I|12fl commutes with scalar poten- 
tials and therefore drops out in Eq. (|ll|l . More generally, 
for arbitrary confining potential, the first term is replaced 
by i(l — V)£ ■ er, resulting nevertheless in no coupling of 
spin to electric fields. The second term in Eq. I|12|l . how- 
ever, allows us to express the coupling of spin to charge 
via the electric field E(t) = (l/e)(V> | VV(r, z, t)\ip ) that 
acts on the quantum dot electron. For the harmonic con- 
fining potential, we obtain 



h(t) 



—e 



(13) 
(14) 
(15) 



The dimensionless field fl(t) describes a combined effect 
of the spin-orbit interaction and electric fields (or more 
generally potential fluctuations) on the qubit. Sl(t) was 
calculated in Ref. for the phonon potential and in 
Ref. |2^ for the shot-noise of a QPC nearby the quantum 
dot. In our case, is merely a classical driving field 
generated by the ac-signal. 

Considering further a constant in space (at least 
on the scale of the quantum dot), alternating elec- 
tric field E(i) — E s'm(uj ac t) of amplitude Eq = 
Eq (cos </>, sin </>, 0) , where </> is the angle of Eq with re- 
spect to the axis x' (see Fig. QJ, we obtain explicitly 
fl(t) = Ho sin(w QC i), with 



—eE , x . x 
S2o = n A_ sm 0, A, 

m e u>Q 



cos 



0,0) 



(16) 



To give an estimate for the amplitude fio in GaAs quan- 
tum dots, we assume A + w A_ w A50 = 8 /im, fkvo = 
lmeV, and E = 10 2 V/cm, which yields Oq ~ 10 -3 . 



The amplitude of the resulting effective magnetic field 
due to EDSR is found from Eq. lfT4T) to be 



5B = 2Bx n . 



(17) 



The maximal amplitude is obtained for B _L J7o, which in 
experiment can easily be arranged for by, e.g., choosing 
B || z. In- plane magnetic fields can also be used, pro- 
vided E(t) is linearly polarized. For example, an elec- 
tric field E(t) aligned with x' generates, according to 
Eq. I|15fl . a dimensionless field fi(f) along y'. In this case, 
B should be chosen along x' for maximal spin-electric 
coupling. 

Using our previous estimate for Qq ~ 10 -3 , we ob- 
tain from Eq. (|17(l that SBq ~ 1 mT for a magnetic field 
B = IT oriented transversely to flo. In principle, the 
dimensionless factor Qq can be increased up to Qq ~ 1. 
However, this requires a specially designed setup, where 
the value of the electron displacement ro in Eq. ap- 
proaches the spin-orbit length Xso- 

Next we remark that fl(t) in Eq. (|15|l can be written 
by the order of magnitude as Q(t) ~ r (t)/Xso- More 
rigorously, we rewrite Eq. I|15|) in the following form 



fi i(*) = E( A so) ii ^(t), 



(18) 



where {X S q) ^ is a tensor of inverse spin-orbit lengths, 



fA-M -I 1/A - 



(19) 



with 1/A± = m e (P ± a)/h and the frame (x',y') was 
used to represent the tensor. For order of magnitude 
estimates, it is useful to introduce the scalar 



1 



A 



so 



(20) 



where ||A S o|| is the Frobenius norm of (A^q)^. In the 
case of Eq. <jT§|) . we have l/Xso — (m e /h)^a 2 + f3 2 . 

Despite the fact that Eq. (|18fl was obtained consider- 
ing the harmonic confining potential as an example, its 
generality suggests that it should remain valid for quan- 
tum dots of arbitrary confinement, provided, to a good 
approximation, the ac-signal merely displaces the quan- 
tum dot parallel to itself by a vector ro(i) as a function 
of time. Note that ro(t) is the only available parameter 
to be compared with Xso in the limit of strong confine- 
ment (Ad — » 0). We extend the class of Hamiltonians 
considered here to any combination of confinement and 
ac-voltage potential that can be rewritten in the form 

U(r) + V(r, t) = U(r- r (t)) + V Q (t), (21) 

where Vo(t) is independent of r. We note that, as be- 
fore, the electron wave function extension A^ is assumed 
to be small compared to the spin-orbit length Xso a t 
each moment in time. Equation l|21|l need not be satis- 
fied exactly. Note that Xd enters only in the definition of 
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ro(t) and does appear alone as a parameter in Eq. (|18fl . 
Therefore, denning r (t) as the average electron position, 
r Q (t) — J r \ip(r, t)\ 2 d 2 r, we expect Eq. (jl%|) to be valid 
to leading order also when the electron probability den- 
sity \ip(r,t)\ 2 changes shape, but the dot size changes 
weakly. 

Equations JTSJ, d}, and (JTHJ form the basis of EDSR 
in quantum dots and can be used to efficiently manip- 
ulate the electron spin by electrical gates. Finally, we 
remark that Eqs. H13fl - I|15|) have been derived under the 
assumption ro <C Xso, an d therefore, can be used only 
for f2o <C 1. In Sec. IVII we discuss the case of f^o ~ 1 in 
more detail. A further assumption in deriving Eqs. 1131) - 
(|15fl was that the frequency spectrum of E(t) lies well 
below the size-quantization energy Hujq. This adiabatic- 
ity constraint is generic to the spin-based quantum com- 
putation^; it guarantees that the electron is not excited 
to higher in energy orbital levels. 



where S 



)n gives the detuning from res- 



onance. The relaxation tensor is diagonal, with 
Txx = Tyy = l/Ta and Tzz = l/Ti, and assumes 
T{ = TijSj . Here, T\ and T 2 are the relaxation and de- 
coherence times in the absence of driving measured in ex- 
periment , 24 i 25 i 31 and S T = -{ng/2\g\)tanh(E z /2k B T) 
is the thermodynamic value of spin, with T being the 
temperature. 

The time-evolution of S in Eq. 1|27|) is simplest in a 
coordinate frame (X' ,Y' , Z'), with Z' |j (6 + u>r), and 
reads 



sin U S 2 



S x ,(t) = Sle-^; 

S Y ,(t) = S°j_e~ t/f2 cos (t^jd 2 + w 2 R + <p \ , 
S Z '(t) = S T + (S° z ,-S T )e- t/fl , (28) 



IV. SPIN DYNAMICS AND COHERENCE 



The electron spin obeys the Bloch equatioi 



.29.30 



(S) = [uz + 6u(t)]x(S)-r(S)+Y, (22) 

where u>z = g/.i B B/h is the Larmor spin-precession fre- 
quency and 6u)(t) = h(t)/H. The spin relaxation tensor 
Tij and the inhomogeneous part Yj are due to the envi- 
ronment and can be derived microscopically 21 ' 29 within 
the Born-Markov approximation. Strictly speaking, Tij 
and Tj in Eq. (|22|l depend also on the driving. In particu- 
lar, Tij acquires, in general, a time-dependent part. How- 
ever, we neglect these effects here since the energy scales 
are well separated. Indeed, from experiments 2 ^ 2 ^! and 
theory, 2 ^ we infer that r 2j ,T 2 - (10 2 - 10 6 )s-\ i.e. 
they are very small, so that the regime Tij , T,; <C 6u> <C 
loz holds. In this regime, the rotating wave approxima- 
tion^ is valid. We consider a completely general driving 
field 

5u(t) = 8u) a sin(u ac t) + Suj b cos(w QC i), (23) 

which can be realized in practice by implementing two 
independent electric fields at the quantum dot site. This 
is, however, by no means necessary for our EDSR scheme. 
The Rabi frequency then reads 

u)r = i {8u) a x n — [<5u>b x n] x n) . (24) 

Here, we assume that u> ac is not far from resonance, i.e. 
\uj ac — uz\ < u ac /2. ln a coordinate frame [X, Y, Z) with 
Z || B, the spin dynamics is approximated as follows 

(S±(t)) « Siftje^"* (25) 
(S z (t)) « S z (t), (26) 

where S± = Sx ± iSy- The spin S(t) obeys a simpler 
(static) Bloch equation 



where Sj_, S z ,, and </> give the initial spin state, (S(0)} = 
S(0) = (Sj_ smcj), Sj_ cos0, S z >), in the coordinate frame 
(X' ,Y' , Z'). Furthermore, the decay times T\ and T 2 
read 



1 



1 



fi S 2 +LU 2 R \T 1 T; 



5 2 , <4 



1 

T 2 



2(S 2 + u% ) V T-j 



2S 2 



To 



(29) 



The stationary value of spin St 
ing order reads 



S(t — > 00) to lead- 



■ tanh {E z /2k B T) . (30) 



S = (d + u R )xS-TS + T, 



(27) 



2\g\ S 2 + (Ti/T 2 )o 

Note that at resonance (S — 0), the right-hand side in 
Eq. (|30|l vanishes. Therefore, in the vicinity of resonance, 
St is determined by the subleading order term, which 
can be obtained from Eq. 130(1 by replacing the numerator 
(S + u>r)5 — > (1/T 2 )[ujji x n}. Measurement of St in the 
presence of driving provides information about the spin 
lifetimes Ti,2- For instance, at resonance the relaxation 
time T\ can be accessed at the leading order of Tij/ujR <C 
1, 

S T (S = 0) = -f^fteaih (E z /2k B T) . (31) 
\g\ 2T x w R 

Finally, we estimate the Rabi frequency ljr using 
Eq. H24|) and the parameters from Sec. IIHI For f2o ~ 
10" 3 , \g\ = 0.44, and B = 10 T we obtain lor - 10 8 s _1 . 
We conclude that, with the present quantum-dot setups, 
EDSR enables one to manipulate the electron spin on 
a time scale of 10 ns, which is considerably shorter than 
the spin lifetimes, for which values between 1 and 150 ms 
(depending on the applied magnetic field) in gated GaAs 
quantum dots have been reported recently^^i 
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V. p 3 -DRESSELHAUS TERMS 

Next we consider the so-called p 3 terms of the Dressel- 
haus spin-orbit interaction^ 3 ., 



HSO = 77 (PyPxPyVx ~ PxPyPx^y) , 



(32) 



where 7 = a c / \/2m\E g is the spin-orbit coupling con- 
stant, with a c (ps 0.07 for GaAs^ 3 -) being a dimensionless 
constant defined in Ref. ^ and E g the band gap. For 
simplicity, we impose here the dipolar approximation for 
the oc-signal, 



V(r,t) = 



■ I dr' ■ E(r',t) re eE(t) 
Jo 



(33) 



Quite remarkably, if the quantum dot potential is har- 
monic, U(r) — u ij r i r ji then the spin does not cou- 
ple to E(t) at the first order of Hso and zeroth order 
of Ez- Indeed, the second term in Eq. vanishes for 
V(r, t) — eE(t) ■ r and S — Hso because of the 
following two identities 

(MLfHso,rMo) = ^o\[L^r,H so }Ho), (34) 
[L^r^Hso] = 0, Vffso(p). (35) 

The latter is specific to Hd in Eq. J^J with a harmonic 
U(r), for which the operator L^r can be expressed via 
the components of p — (e/c)B z x r. Note that, generally, 
[p, Hso] = (e/ c) [B z xr, H S o] for any H S o that is a func- 
tion of only p = (px,Py)- Thus, for a harmonic confining 
potential, one is left with the same dominant mechanism 
as considered above for the "linear in p" terms. Ex- 
panding in terms of the Zeeman interaction, we recover 
Eqs. ifDfy and ifTljl with fl(t) given now by 



fli(t) 



m e ojQ 



(36) 
(37) 



where (Xg Q )ij is a tensor of inverse spin-orbit lengths, 
and as before we consider U(r) = m e uj 2 r 2 /2. For Hso 
in Eq. ((321) . we obtain explicitly (Xg )ij = k'yuo m e.^ij 
and Qi(t) = —jem e Ei(t) /4u>o- To estimate the strength 
of the resulting EDSR, we note that 7 ~ (3d 2 /ft 2 , and 
therefore the amplitude of hit) = 2gfj,sB x tl(t) is down 
now by a factor d 2 /X 2 <C 1 compared to the p-terms. 

Next we consider a quantum dot with non-harmonic 
potential U{r) and show that the p 3 -terms in Eq. il-il'l) 
give rise to a spin-electric coupling proportional to the 
cyclotron frequency uj c — eB z /m e c. Since ftio c differs 
parametrically from Ez (Ez/ftuJ c = gm e B/2mB z ), the 
p 3 -terms can be as significant as the p-terms, provided 
Ez/ftuJc < d 2 /X 2 , which is realistic for GaAs quantum 
dots. Note that for the p-terms no spin-electric coupling 
proportional to ui c arises at the first order of Hso- We 



thus leave out Lz in Eq. 1(10(1 and consider a confining 
potential U(r) that differs from a harmonic one by a 
function W(r), 

U{r) = Uijnrj + W(r) = U H (r) + W(r), (38) 



where real coefficients and W(r) — 0(r 3 ). While 

in general W(r) need not be small compared to Hh = 
p 2 /2m e + Uh{t), in the following we expand the denom- 
inator of Eq. (|10|l in terms of W <C Hjj, considering 
therefore only small deviations of U(r) from harmonic 
potentials. Then, using Eqs. 113(1 . and (|33|l . we 

obtain at leading order in lu c 



dui r da,, 



(39) 



.=0 



where we set uj c — > in the right-hand side of Eq. 1(39(1 
after evaluating dSu/doJc in the symmetric gauge, with 
Sh defined as [Hh,Sh] — Hso- The linear relation- 
ship between hi and oj c holds for lu c <C luq, where 
ujo = 2y/det(u) /m e Ti(u) . In Eq. (|3*9"|) . the perturbation 
W(r) enters via the function R(r) defined as follows 



Ri(r) 



i 



dW{r) 
'v dn ' 



(40) 



Note that (R/r) ~ WoXd/fkJoXw is the small parameter 
of our expansion in terms of W(r), with Wq and Xw < 
Xd being, respectively, the characteristic amplitude and 
length scale of the variation of W(r) over the quantum 
dot size. It is important to note that the antisymmetric 
part of W(r) drops out in Eq. because Hso is also 
antisymmetric with respect to r — ► — r. 

Next, as an example, we consider £/jj(r) = m e ojQr 2 /2 
and W(r) = rjr 4 , and obtain 



\h{t) ■ a 



ejrjh 2 uj c 
9m e WQ 



(E y (t)a x + E x (t)<Jy). 



(41) 



Here, we have used the deformation quantization the- 
ory^ which allowed us to considerably simplify the 
derivation of Eq. ((41(1 by performing most of the calcu- 
lation in classical mechanics and only at the final stage 
come back to quantum mechanics. We have also carried 
out a fully quantum derivation of Eq. 1(41(1 and recovered 
the same result. 

To estimate the strength of the resulting EDSR, we 
note that h ~ ftw c (Xd/Xso)(eXdEo/ftu!o){R/r), where 
X so = (« A|/[0.01nm] for GaAs) is the 

spin-orbit length of the p 3 -terms and the parameter 
{R/r) ~ W(Xd)/ftoJo characterizes the deviation of the 
quantum dot confinement from harmonic. In practice, 
{R/r) can be as large as unity, but here we assume 
{R/r) ~ rjX d /fto-!o = 0.1. For an electric field with am- 
plitude Eq = 10 2 V/cm and a GaAs quantum dot with 
ftujQ = 1 meV, we obtain the equivalent of an ac mag- 
netic field 8B(t) = h(t)/gfiB that has an amplitude 
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SB r* lmT at B z = IT and |g| = 0.44. In contrast 
to the previous mechanism, 8B(t) can have here also 
a finite longitudinal component 8B\\(t) = n(n ■ SB(t)), 
which however vanishes if B \\ z. 

Finally, we note that the p 3 -terms can also be relevant 
for spin relaxation in quantum dots with non-harmonic 
confining potential. Of course, the magnetic field has 
to have an out-of-plane component for this spin-electric 
coupling to dominate over the one considered in Sec. IIIII 

VI. DISCUSSIONS 

The coupling of spin to electric fields that we have de- 
rived above can be used in a variety of ways to access and 
manipulate the electron spin in experiments. The effec- 
tive Hamiltonian in Eq. has the same form as the 
Hamiltonian of an ESR effect. This shows that ESR and 
EDSR arc mutually interchangeable and the choice of the 
effect to be used depends on the particular experimental 
setup. In GaAs quantum dots, the spin-orbit interac- 
tion is weak enough to ensure long coherence times and, 
at the same time, strong enough to allow room for spin 
manipulation on an experimentally accessible timescale 
of ~ 10 ns. Much shorter timescales can be achieved 
in In As quantum dots, because of a much stronger spin- 
orbit coupling and a larger electron g-i actor. In con- 
trast, the EDSR effect is of little use in materials with 
very weak or nearly absent spin-orbit interaction, such 
as, e.g., carbon- nanotube quantum dots. To make order- 
of-magnitude estimates easier, we draw an analogy be- 
tween the EDSR and ESR effects in terms of the partic- 
ular way the B and E fields couple to the electron spin 
S= (h/2)(T. 

We recall that the ESR effect occurs as a result of 
the Zeeman interaction of the electron spin with an ac 
magnetic field. It is convenient to write this interaction 
in the form of a magnetic dipole interaction, 

H ESR = -n ■ B(t), (42) 

where B(t) is the ac magnetic field and fi is the electron 
magnetic moment, 

1 . , 

V = --^g^Bcr, (43) 

where g is, in general, a tensor, see Eq. (|A18I) . 

By analogy with the ESR effect, the EDSR effect can 
be viewed as arising from an interaction between the ac 
electric field E(t) and a spin-electric moment v. The 
respective spin-electric interaction is then analogous to 
Eq. and reads 

Hedsr = -v ■ E(t), (44) 

where the spin-electric moment u is due to an inter- 
play between the spin-orbit interaction and some time- 
reversal breaking interaction, such as the Zeeman inter- 
action. This analogy is not complete. Equation (|44|l is 



valid only for ac electric fields E(t) that oscillate around 
zero, whereas Eq. (14211 holds also for static B-fields. The 
reason why a static electric field E cannot be used in 
Eq. (|44|l will become clear after we explain the origin of 
v in Eq. (I4U) . 

The spin-electric moment v arises because the dipolar 
transitions in the quantum dot become allowed, e.g., for 
the ground state 

MnHVcu) y^O. (45) 

The electron charge density operator p{r) = —eS(r — 
r e i), where r e i is the electron coordinate, acquires spin- 
dependent terms in the transformed basis, 

|^ s ) = e- s |^ n )| Xs ), (46) 
p(r) = e s p(r)e- s , (47) 

where e~ s is the transformation used in Section lllll and 
studied in detail in Appendix [S] One can present p(r) 
as a sum of two terms, 

p(r) = p(r) + 6p(r), (48) 

where p(r) is spin independent and Sp(r) is proportional 
to the spin. Then the spin-electric moment can be writ- 
ten as follows, 

v = J rSp(r)dv, (49) 

where dv is the elementary volume of integration. Equa- 
tion l|49l) unveils the physical meaning of the spin-electric 
moment v: due to the mixed spin and orbit nature of the 
electron density, the electron spin couples to the first mo- 
ment (dipole moment) of the electron. 

While oscillating around an equilibrium position, the 
electron produces a time-dependent dipole moment, part 
of which is proportional to the electron spin. Obviously, 
in a static electric field, one can set to zero the electron 
dipole moment, because the new electron position can 
be taken as the equilibrium one. Therefore, only the 
change of the moment as a function of time has a physical 
meaning for single electrons (as for any other monopoles) . 
In contrast, the electron magnetic moment fi couples to 
static magnetic fields, because one can view fi as arising 
from a pair of Dirac monopoles of opposite signs, for 
which the relative distance between them has an absolute 
meaning. 

Equation 1|49|) is written in a very general (operator) 
form. After taking the expectation value in the orbital 
ground state IV'o); we obtain 

v = -Ua, (50) 
d 

Pij = 2e—{ipo\e s ne~ s \ijj ), (51) 

where the derivative with respect to Uj is defined as a 
usual derivative of an expression that is linear in cr. Note 
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that Eq. (|50(l is analogous to Eq. I|43[) where the role of 
gfiB is played by the tensor v. Using Eqs. (|T3)l - lfH)l) we 

obtain for the linear in momentum spin-orbit interaction, 



(52) 



Similarly, for the p 3 Dresselhaus terms we obtain from 
Eq. gTJ) 



9m e Wg 



1 

1 



(53) 



where we use the coordinate frame (a;, y) to represent 
the tensor. Note that in both cases Vij is proportional 
to the magnetic field (or one of its components). The 
spin-orbit interaction produces no spin-electric coupling 
at B = 0, because of the time-reversal symmetry of spin- 
orbit interaction. 

The analogy between fj, and v is also seen in the pair- 
wise interaction between spins in separate (not tunnel 
coupled) quantum dots. 34 For an unscreened Coulomb 
interaction between electrons, the spin-spin interaction 
is analogous to the magnetic dipole-dipole interaction^ 

„ Ui ■ i^rf - - 3(i/ 4 ■ Tij){v s ■ ry) 

ff dd = 2^ " ^5 ' ( 54 ) 

where r,j = r, — Tj is the distance between two quantum 
dots (yij S> Xd) and k is the electric permittivity of the 
material. For further detail and a microscopic derivation 
of Eq. (|54|l we refer the reader to Ref. 0- 

Next we discuss the limitations of our theory. 
Throughout the paper, we have assumed that the spin 
orbit interaction is weak compared to the dot level spac- 
ing, or, in other words, that Xd/Xso "C 1. This assump- 
tion allowed us to use the perturbation theory to find 
the unitary transformation exp(— S), see AppendixEl Of 
course, this restriction was not necessary, if, e.g., we were 
to apply numerical methods to diagonalize the Hamil- 
tonian. In particular, note that Eq. Ij49(l is meaning- 
ful whenever the unitary transformation exp(— S) exists. 
The latter is always the case, including also extended 
states. Our perturbative results are qualitatively correct 
for Xd < Xso and can be used in experiments to estimate 
the strength of the EDSR effect. The case ^> Xso is 
more seldom and requires a separate theoretical investi- 
gation. 

As a second limitation, we would like to mention the 
adiabaticity criterion. In Sections IlIII and Ivl we have de- 
rived effective Hamiltonians for the low energy subspace 
of the quantum dot Hilbert space. For the validity of this 
effective description, it is important that the switching 
(on and off) of the effective interaction occurs on a time 
scale that is larger than the inverse level spacing in the 
quantum dot. Obviously, this criterion excludes applica- 
bility of our theory to extended states. In practice, how- 
ever, the finite temperature T imposes a more stringent 
criterion on the confinement energy, h 2 /m e X d 3> ksT. 



The third limitation of our theory is a small amplitude 
of oscillation of the quantum dot, r /Xso "C 1. We have 
shown in Section lTTTl that the EDSR effect is proportional 
to this small parameter. Thus, for breaking the time- 
reversal symmetry by the Zeeman interaction we have 
(by order of magnitude) 



U)R ~ ui z 



TO 

Xso ' 



(55) 



where ojz = Ez/h. Similarly, for breaking the time- 
reversal symmetry by the orbital B-field effect (Sec- 
tion 01 , we have 



Xso 



(R/r), 



(56) 



where {R/r) is the small parameter of deviation of 
the quantum dot confinement from harmonic. We re- 
mark that our theory remains qualitatively valid also for 
r o/Xso ~ 1 an d for (R/r) ~ 1. Beyond these limits, we 
do not expect the Rabi frequency to grow indefinitely 
The Rabi frequency is bound in the case of Eq. (|55l) by 
< cjz, and in the case of Eq. I|5t)|l by lor < oj c . We 
conclude that, by designing quantum dot setups that al- 
low for large oscillation amplitudes ro < Xso, the EDSR 
effect can be strongly enhanced, beyond the numeric es- 
timates made in Sections Mil IIVI and ^ 

In conclusion, the EDSR mechanisms presented above 
provides a means of implementing local electrical control 
of electron spins in quantum dots. 
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APPENDIX A: SCHRIEFFER- WOLFF 
TRANSFORMATION AND FINE STRUCTURE 

In this Appendix, we first work out the Schrieffer- 
Wolff transformation to the third order of perturbation 
theory and for a general weak perturbation. Then, we 
consider an example Hamiltonian and use the Schrieffer- 
Wolff transformation to partly diagonalize the Hamilto- 
nian. Finally, we analyze the fine structure of the trans- 
formed Hamiltonian and complete its diagonalization by 
an additional unitary transformation. 

As in standard perturbation theory, we consider a 
Hamiltonian H = Hq + Hi, where Hi is a weak per- 
turbation with respect to Hq. For the matrix elements 
of Hi , we assume 



(n\Hi\m) = 0, for E n 



E„ 



(Al) 



(n\Hi\m)<$:E n -E m , for E n + E m , (A2) 
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where |n) and E n are, respectively, the eigenstates and 
eigenvalues of H , and are obtained from H \n) = E n \n). 
The projector V , defined as follows 

VA= J2 A nm\n){m\, VA (A3) 

projects onto the diagonal or degenerate part of Hq. In 
the particular case, when the spectrum of Hq is non- 
degenerate, V assumes PA = J2 n A nn \n)(n\, WA. From 
Eq. 1A1(> and the definition (|A3f) , it follows that 

VH 1 = 0, (A4) 
VH = H . (A5) 

Next we look for a unitary transformation that brings 
the Hamiltonian H = Hq + Hi to a partly diagonal form, 

H = e s ' {H + Hi ) e- s ' = H + AH, ( A6) 

where the operator AH obeys VAH = AH. Here, 
S 1 = —S'' is the transformation matrix. The unitary 
transformation in Eq. (|A6|) is called the Schrieffer- Wolff 
transformation.— We expand <S" and AH in terms of the 
perturbation H\\ 

S' = + S' {2) + S' (3) + . . . , (A7) 
AH = AH^ + AH^ + AH^ + ..., (A8) 

where the superscripts give the order of perturbation the- 
ory. Substituting Eqs. l|ATjl and (|A"8|) into Eq. <jA"6l) . we 

find a set of equations for S' , 

[H ,S'W] = H U (A9) 
[H ,S'W] = f ffr], (A10) 

[H ,S'W] = ^[S^\H 1 ] + ^[S'^,[S'^,H 1 }} + 

Q[S'W,V[S'W,H l ]], (All) 

where Q = 1 — V . It is important to note that S is 
defined in Eqs. (jA^|) - (jAll|) up to terms VM, where M 
is arbitrary. Such terms drop out on the left-hand side 
in Eqs. ifXflTl-lfXTTTl because [H ,VS r \ = 0. Thus, VS' 
can be chosen arbitrarily, which shows that there are 
infinitely many transformation matrices S' that satisfy 
Eq. I|A6(I . For simplicity, we choose VS' = and address 
the fine structure of H = H + AH later on. For the 
operator AH , we obtain 

AffW = 0, (A12) 

AH (2) = ^[S'W.ffi], (A13) 

AH^ = (A14) 

Introducing the Liouvillean L n : L a A — [H ,A], VA, 
we can formally solve Eqs. l|A9p - (|All|) one by one. For 



example, the transformation matrix at the lowest order 
reads S^ 1 -* = QLq Y H\. For AH, we recover then the 
perturbation theory expansion in a more familiar form, 

AH = -VHiL^Hx + VHxt^Hxt^Hx + ..., (A15) 

with the usual convention, VL^ 1 A = 0, \/A, adopted. 

Next, we remark that the fine structure of H = 
H + AH can be addressed in each particular case by 
means of degenerate perturbation theory. As an exam- 
ple, we consider here the Hamiltonian H = Hq + Hi, 
with Hq = Hd + Hz and Hi = Hso ■ Here, Hd is given 
in Eq. (|5j), with U(r) = m e WQr 2 /2, Hz is given in Eq. Jfjjl, 
and Hso is given in Eq. Q. Using the transformation 
matrix S' = S, with S given in Eq. (|12f> . we obtain a 
diagonal Hamiltonian, H — Hd + Hz, at the first order 
of Hso- At the second order of Hso, however, a fine 
structure in the energy spectrum arises. At B = 0, the 
transformed Hamiltonian reads 

where £ z = — i(xd/dy — yd/dx) is the electron rotational 
momentum and Aso = 2m e (/3 2 — a 2 ). The Kramers 
doublets are identified, in this the pairs of states 

with quantum numbers {l z , a z ) and (—£ z , —cr z ). For l z > 
0, the two-fold orbital degeneracy is lifted and a splitting 
Aso^z arises. Note that the ground orbital state, which 
has l z = 0, remains doubly degenerate in this case. 

At B 7^ 0, the fine-structure interaction in Eq. (|A16|) 
is modified by both the Zeeman energy Ez and the 
cyclotron frequency uj c . For simplicity, we omit terms 
~ AsoEz /h/jjQ, but keep terms ~ Aso^c/^o, assuming 
that Ez <C hjj c . Then, the Hamiltonian (|A16|I acquires 
two extra terms 

EZ . ASO rpy 2 

— n • o- + —^-cr z Vr , (A17) 

where A = \Jh/m e Lo c is the magnetic length and we use 
the symmetric gauge, A(r) = B z {—y/2, x/2, 0). The last 
term in Eq. I|A17(I can be viewed as a renormalization of 
the electron g-factor. Allowing for an anisotropic Zeeman 
interaction, 

ffI ff = 5^E^ B ^ ( A18 ) 

ij 

we obtain that the tensor gij is diagonal in the main 
crystallographic frame, with 

9xx = 9yy = .9i 

g zz =9 + ^^(Ur 2 \^n), (A19) 

where m is the electron mass in vacuum and ip n is the 
electron orbital state. For the ground orbital state, the 
corrected (/-factor reads g zz = g + mAsol m J^>, where 
u> = ^/lUq + ojI/A. Note that the sign of the correction is 
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given by the sign of 1 — a 2 contained in Ago- The spin 
quantization axis does not, in general, coincide with the 
magnetic field direction n and is given by the following 
unit vector 



n + (n z 



VT+C(2 + CK : 



(A20) 



where £ = (g zz — g)/g. An additional unitary transfor- 
mation can be used to diagonalize the 2x2 blocks of 
Zeeman-split Kramers doublets, 



Hf 



2 s " Hf 



-E z n ■ <x 



(A21) 



where E z = E Zy /l + C(2 + O n l is 

the renormalized Zee- 
man energy and 



-s" 



n ■ n . £[n x n z ] ■ rr 



1 — n ■ fi 



So far, we have considered a given orbital state ip n , for 
which the tensor is given in Eq. I|A19|I . The transfor- 
mation above is also valid in general, provided £ is un- 
derstood as a diagonal operator, ( — (mAso/ g^'Pf 2 ■ 

We summarize by mentioning that the unitary trans- 
formation in Eq. i|A6|l can, in principle, be adjusted to 
give a fully diagonal H = Hq + AH, i.e. we had not to 
require VS' — in the first place. However, in practice, 
it is more convenient first to apply the non-degenerate 
perturbation theory, Eqs. (jA9jl - (|A15p . and then, at the 
end, complete the diagonalization of Hq + AH by a sec- 
ond unitary transformation. The latter is specific to each 
particular case and amounts, in general, to applying the 
degenerate perturbation theory. For the sake of simplic- 
ity, we shall refer to S in the main text of the paper as 
to the full transformation matrix, despite the fact that 
the respective unitary transformation comes, in practice, 
as a product of two unitary transformations. Thus, we 
denote the product e~ s e~ s by e~ s in the main text. 



(A22) Finally, we remark that e~ s " ss 1 + 0{H 2 S0 ) 
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